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Abstract 

We extract the matrix elements of four quark operators Ol s relevant to the B s and B s 
life time difference from QCD sum rules. We find the vacuum saturation approximation 
works reasonably well, i.e., within 10%. We discuss the implications of our results and 
compare them with the recent lattice QCD determination. 
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1 Introduction 



The recent results on CP violation in Bd - Bd mixing have been reported by the BaBar and Belle 
Collaborations |l| in the ICHEP2000 Conference. More experiments on B physics have been 
planned in the present and future B factories 0. Theoretical efforts to improve predictions 
and reduce uncertainties are expected and needed. It is well-known that mixing in neutral B 
meson systems provides a good place to examine CP violation as well as flavor physics in the 
standard model and beyond. For example, the mass difference between the mass eigenstates 
of neutral Bd meson, AM# d , gives an important constraint on CKM matrix element V t d and 
the first indication of large mass of top quark. Similarly, the mass difference between the 
mass eigenstates of neutral B s meson, AMb s , which will be precisely measured in the near 
future would give an valuable constraint on CKM matrix element Vt s . The another important 
observable for mixing in neutral B meson systems is the lifetime difference between the mass 
eigenstates of neutral B mesons, AT Bd or AT Bg . The ratio |V ts /V^| 2 can be extracted from 
the measurement of ATb s ||. The width difference of Bd mesons is CKM suppressed and 
consequently not easy to be observed. In contrast, for B s mesons the width difference is large 
enough to be measured H and has been recently measured || with low statistics. Hopefully, 
it will be measured with high statistics in the near future. 

As usual, The light B^ and heavy mass eigenstates are defined by 

\BL> H )=p\B° s )± q \B° s }, 

where \B®) and \B®) are the flavor eigenstates. The mass difference and the width difference 
between the physical states are given by 

Am = Mh — M L , AT = T H - T L . 

Because | r i2 1 -C |Mi 2 | for B s mesons ||, to the leading order in ^12/^12!, Arris = 2|Mi 2 |, ATb = 
23?(Mi2Pi 2 )/|-^i2|- Neglecting very small CP violating corrections, the width difference for B s 
mesons in SM has been given || [7| 

(y) fl = ( ai^Mev ) [°- 006fi K) + °- 150 ^K)- 0.063], (1) 

where fg a is the decay constant of B s , B and Bs are the bag parameters related to the four 
quark operators Ol and Os (see below). These hadronic quantities need to be calculated by 
non perturbative methods such as lattice, QCD sum rules, Bethe-Salpeter approach, etc. 

The similar quantities related to B® — B^ mixing have been estimated by Narison et al within 
the traditional QCD sum rules approach [|J at a s order. Their conclusion is that the vacuum 
saturation values Bb — Bb* — 1 are satisfied within 15%. Their sum rules are constructed 
through two-point correlation functions and depend on some phenomenological assumptions. 
In this letter we shall calculate the matrix elements of four-quark operators relevant to the B s 
meson lifetime difference through QCD sum rules in HQET. The sum rules are constructed with 
three-point correlation functions. Our calculation is carried out at the leading order in 1/rrib 
expansion in HQET for simplicity. In ref.|| the effects of condensates are absorbed like other 
factorizable corrections into the contribution to fs and the available result of H (though not 
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explictily said) has been used as the effects are small compared to the perturbative corrections. 
In our sum rules the nonperturbative contributions of condensates are explicitily included and 
the numerical results confirm the smallness of these corrections (see below). 



2 Theoretical formalism 

We employ the following three-point Green's function, 

r>V) = z 2 f ^e* fc '^(0|T^ , (2) 



where lo = v • k, u' = v • k'\ is the b-quark field in the HQET with velocity v. And Ol,s 
denotes the color-singlet four quark operators. They are 

L = 6 7 M (l-7 5 )sMl-75)s, (3) 

O s = 6(1 - 75 )s6(l - 76 )s, (4) 
In terms of the hadronic expression, the correlator in Eq. (|2|) reads 

r ot i\ F b s (B S \Q LS \B S ) 

r ( W , u, ) = — (A _^ (A _^ + resonances , (5) 

where A = itlb — mi, and Fb„ is the B s decay constant in the leading order of heavy quark 
expansion defined as [] 

(0|s(0)75^ } (0)|^) = -i^F Bs (7) 

In order to eliminate the contribution from the non-diagonal single pole terms and suppress 
the continuum contribution in Eq. (5), we make double Borel transformation to the correlator. 
The transformation is defined as 

A= Hm l lm t^(jXt^L(AY. (8) 

-c^oo -uj'^oo n - \ duJ J m] \ duJ 'J 

n — > oo m —y oo 

f = fixed f' = =^- fixed 

n m 

There are two Borel parameters f and f, which appear symmetrically, so f = f — 2T are 
taken in the following analysis. 

On the other hand the correlator can be calculated at the quark gluon level. For example 
for Ol we may rewrite the right hand side of Eq. (fj) as 



"Note that Jb s in Eq. (1) is denned by 

(0 | a r 75 b \B Q S ) = -if BsP » (6) 
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- 2 J dxdye M - x+lk y{-Tr[ l5 • zSf(x) • - 75) • ^(-y) • TirfG/) • 7^(1 - 7s) • «^ TO (-a:)] 

+Tr[iSr(-x)] ■ 75 ■ AS^s) • Y(l - 7 5 )]Tr[z5f(-y)] ■ 75 ■ i5?(») ■ 7^(1 - 7s)]} (9) 

where iS J s n (x) is the full strange quark propagator with both perturbative term and condensates, 
i,j etc is the color index. iS™ J '(x) is the leading order heavy quark propagator which has very 
simple form in coordinate space: 

poo 

iSt J (x) = 5 ij / dtS(x-vt) (10) 
Jo 

Note the structure of color flow is quite different for the two terms in Eq. ([jp. For the 
perturbative part the first and second term is proportional to N c and N%, respectively, where 
N c = 3 is the QCD color number. In the limit of N c — > 00, the second term dominates! As 
shown below, the non-factorizable contribution in Fig. 1 d, f and g has different color structure 
from the factorizable terms in Fig. la, b, c and e. The condensates up to dimension six are kept 
in our calculation. We also expand the strange quark propagator and keep perturbative term 
of order 0(m s ). The calculation is standard and we simply present final results after making 
the double Borel transformation. 



3 Duality Assumption 

We may write the dispersion relation for the three-point correlator T(u,uj') as 

„. , s 1 r°° , r°° , , lmU(u,i/) , , 

II(uj,u') = — dv dv 1 - ± ' ; - . 11 

v ' ; TV 2 Jo Jo {v - u){v' - uj') v ; 

In order to subtract the continuum contribution, we have to invoke quark hadron duality 
assumption and approximate the continuum by the integral over the perturbative spectral 
density above a certain energy threshold ui c . 
With the redefinition of the integral variables 

v + v' 



2 ' 
1 



v — v 



(12) 



the integration becomes 



dv \ dv... = 2 dv + dv (13) 

Jo Jo 



It is in v + that the quark-hadron duality is assumed JlO], |TT], [12] 



1 * 

, , 2 f°° 1 f v+ 1 ImIT(z/, v'\ , 
higher states = — - / dv + / dv--. r-. r . (14) 

7T Jlu c J-u + \V — UJ){V' — U)') 
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This kind of assumption was suggested in calculating the Isgur-Wise function in Ref. jlT] and 
was argued for in Ref. [f[^| . As pointed out in |T0| , [12| , in calculating three-point functions 
the duality is valid after integrating the spectral density over the "off-diagonal" variable z/_ = 



\[y — v'\ Such a duality assumption is favored over the naive one: 



higher states = — / du dv '- 

Ti Ju) c Ju) c 

4 QCD sum rules 



ImIT(z/, v'\ 



v 



(15) 



The spectral density Pl,s( s i, $2) of the perturbative term reads 

iV c (iV c + 1) 



PL{S 1 ,S 2 ) 



PS{Sl,S 2 ) 



2vr 4 
iV c (2iV c - 



4vr 4 



-s 1 s 2 [sis 2 + m s (si + s 2 )\ 



-sis 2 [sis 2 + m s (s 1 + s 2 )} 



(16) 
(17) 



The sum rule for (B s after the inclusion of the condensates and the integration 
with the variable z/_ is 



^f-(B.\0 L \B a )exp (~ 



N C (N C + 1) /■<*= 



7T' 



{ / (ii/exp 



V 

T 



.16 

— H — mj' 
15 3 



+ 3^(1 
A^ 2 - 1 , 



rrir 



[2T 2 (^G 2 ) + a s m^] 



256tt 4 

where a s = — (27r) 2 (ss) and we have used the factorization assumption for the four-quark 
condensates. Similarly we have 

N C (2N C 



(B s \O s \B s )exp 



T. 



2tt 4 
A 2 - 1 



.16 k 8 4 , 
15 3 J 



rrin 



, . + -m s a s T 2 + — } 
64T 2 ' 6 288 J 



512tH 



[2T>(giG>) + a s mlT] 



(19) 



We want to emphasize that in Eqs. (|T8|), ([19]) the terms with color factor A C (A C + 1) and 
N C (2N C — 1) come from the factorizable diagrams in Fig. 1 a, b, c and e. The non-factorizable 
contribution has a color factor N \ 1 which comes from the summation over color factor, 
Tr[4r4^-] = N \ l , in Fig. 1 d, f and g. A second observation is that the factorizable terms 
are all positive while nonfactorizable pieces are negative. 

Now we turn to the numerical analysis. The decay constant and binding energy of the B s 
meson at the leading order of heavy quark expansion can be obtained from the mass sum rule 



<b s exp 







(-4) 


= 3 1 




8tt 2 Jo 



/ dss(s + 2m s )e~ 
Jo 



-s/M _ 



< SS > (1 - 



AM 2 ' 



(20) 



Note M = 2T, so = 2u c . We have not included a s corrections in Eq. because they are 

also neglected in the sum rule for (B s \Ol s\Bs) (HD"(H)- The values of the parameters are 
calculated to be F Bs = (0.49 ± 0.1) GeV 3 / 2 , A = (0.68 ± 0.1)GeV with the threshold s to be 



(2.2 ± 0.3)GeV and the Borel parameter M in the window (0.65 — 1.05) GeV [13]. Numerically 
we use the following values of the condensates, 

(ss) ~ -0.8 x (0.23 GeV) 3 , 

(g 2 G 2 ) ~ 0.48 GeV 4 , (21) 

(gscr^G^s) = m 2 (ss) , m 2 ~ 0.8 GeV 2 . 

For the strange quark mass we use m s = 0.15 GeV. 

In order to minimize the dependence of the parameters we divide Eqs. (0, [L9|) by Eq. (|20|) 
to extract the matrix elements, the variation of which with uj c and T are given in Fig. 2 and 3. 
The sum rule window is T = (0.2 — 0.5) GeV, which is almost the same as that in the two-point 
correlator sum rule. We obtain 

(B s \O l \B s ) = (0.85 ± 0.20)GeV 4 , (22) 

\(B s \O s \B s }\ = (0.55 ± 0.15)GeV 4 (23) 

where the central value corresponds to T = 0.3GeV and u c = l.lGeV. The uncertainty includes 
the variation with T and uo c . The bag parameters B and B$ are defined by 

8 5 M 2 

(B s \0 L \B s ) = -f 2 M 2 B B, (B s \O s \B s ) = --f Ba M Bg ^ + \ )2 B S . (24) 

and they can be directly obtained from Eqs. (22), (23) and (24). 

The ratio of these two matrix elements is very interesting. We divide Eq. (19) by Eq. (18) 
to extract the numerical value of the ratio. In such a way the dependence on the the Borel 
parameter and the continuum threshold is minimized as can be clearly seen in Fig. 4. Within 
the accuracy of QSR the curve in Fig. 4 is flat. The ratio is practically the same in the working 
region of T and u c . It reads 

\(B s \O s \B s )\ , . 
11 = (Bs\Ol\Bs) = (0 ' 63 ± °' 13) (25) 

In our numerical calculation, the contribution of the perturbative term is about 45 — 65% of 
the total contributions in preferred Borel variable region. We have used factorization approx- 
imation for the four quark condensates in numerical calculations. This may introduce some 
uncertainty. We may introduce a scale factor k to indicate the deviation from the factorization 



approximation as in [14]. In our calculations of sum rules the 1/M& corrections in HQET have 
not been included in , which may bring a deviation from the numerical results of the matrix 
elements. However, for the ratio of the two matrix elements, we expect little change to the 
above analysis. Our numerical results are not sensitive to the mass of the strange quark. Ac- 
tually, the effects due to the strange quark are very small so that the results for B s are almost 
the same as those for B&. 
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We now give a remark on the usual factorization assumption. In our Feynman diagram 
(Fig. 1) calculations, the contributions of nonfactorizable diagrams are around —6%, —7% for 
(B s \Os\B s ) and (B s \Ol\B s ) respectively, which means that the factorization approach works 
well even though our calculations are limited to the leading order in the l/m& expansion in 
HQET. That is, the conclusion in Ref.|| remains unchanged when the nonperturbative con- 
densate contributions are taken into account. If one considers a s corrections, there is only 
one nonfactorizable perturbative diagram, in which the gluon line in Fig. 1 f is connected, 
in the fixed-point gauge in the leading of 1/M& expansion. However, radiative corrections are 
generally of high order 24 compared to the leading order, and the fact is that the perturbative 
term is about 45 — 65% to the whole contribution, so the contribution from the diagram can 
be neglected compared to those in Fig. 1 (d), (f) and (g). The case here is different with that 
in the calculation of matrix elements of four-quark operators, relevant to the life time differ- 
ence between heavy mesons, where the flavor changes AF = 0. In that case, the perturbative 
contribution vanishes [15|1, and we can't predict naively how large the radiative correction is 
compared to the nonperturbative terms. 



5 The B s and B s decay width difference 

The complete expression for AT Ba with short- distance coefficients at NLO in QCD is given by 



lQn 2 B(B s ^Xeis) f 2 B M B 2 



g(z)vQCD 



ml 



(26) 



Ml 



G{z)-B + G s {z) 



B s + VT^4z <5i 



/m 



where 



G(z) = F(z) + P(z) and G s (z) = -(F s (z) + P s (z)). (27) 

and F, P, F$, Ps can be found in Ref. @. We eliminated the total decay rate T Bs in favor of 
the semileptonic branching ratio B(B S — > Xeu), as done in ||. This cancels the dependence of 
(AT /T) on V c b and introduces the phase space function 



g(z) = l-i 

as well as the QCD correction factor [[16] 



\z + 8^ 3 - z 4 



12z 2 In*, 



Vqcd = 1 - 



2a s (m b ) 
3tt 



31 



71 



(2f 



(29) 



One can also express the width difference as 



Ar 



r B Am Bi 



™>B S 



(exp.) T/ 2 
Vis 



V, 



td 



K- [G(z)-G s (z)TZ(m b )+5 1/m ) £ 



(30) 
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where 



K is the Eq. (7) in ref. JT7 



fB d \fB Bd 



$l/m = fB 3 M B 3 $l/m, 



(31) 



(32) 



and the 5i/ m represents the 1 /m b corrections and can be found in Ref. 0. 

It is clear from the above equation that besides the ratio 1Z of the matrix elements of four 
quark operators, which are those we have calculated in the paper, we only use the experimental 
-B^-meson mass difference, which is known with a tiny error [ITS 



(Am Bd f^ 



0.484(15) ps 



(33) 



and another ratio of hadronic matrix elements, £, which is rather accurately determined in 



lattice simulations 19, 20 



As it is well known, the quantities in Eq. |I] are calculated at the scale 0(m b ), while our 
result Eq. (|25"D is calculated at the hadron scale [ihad- Therefore, We have to consider the 
renormalization scale dependence of those four-quark operators. The anomalous dimension 
matrix of these operators has been given in Ref. J7|. Using the anomalous dimension matrix 
and following the standard way, we obtain the scale dependence of 1Z 



n{m b ) = imn(fx had ) + om. 



(34) 



where TZ^had) is defined by Eq. (|25f) . To obtain the numerical result, m b = 4.8 GeV and 
fJ'had = 1-0 GeV have been used. It is obvious from Eq. ( |3"4| ) that the result depends on the 
renormalization scale heavily. 
Numerically we have 



Ar, 



[(0.5 ± 0.1) + (13.8 ± 2.8)n(m b ) + (15.7 ± 2.8)(-0.55 ± 0.17)] x 10~ 2 

= (7.0 ±0.8) x 10~ 2 



(35) 



Clearly such a life time difference is compatible with existed literatures. It's interesting to 



compare our result to the two recent lattice QCD calculation: 



Ar, 



10" 2 in Ref. 21 and 



(4.7 ± 1.5 ± 1.6) x 10" 2 in Ref. |7 



= (10.7±2.6±1.4±1.7) x 
In Eq. (j35|) the numerical 



value of Si/ m , which corresponds to the l/m b correction in the short distance expansion of the 
operator product H e ff(x)H e ff(0) [^], has been taken as -0.55|17]]. If it is taken as -0.30, one has 
Ar/r = 10.9 x 10" 2 , larger than 7.0 x 10" 2 , while in the case of Ref. [|T7|, AT/T would remain 



in the 10% range with the change from -0.55 to -0.30. That is, the sensitivity to the final term 
in Eq. (^), i.e., the l/m b correction, increases in our result. Without a good control of the 
correction, a precise determination of the lifetime difference is impossible. 
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6 Conclusion and discussion 



In summary we have calculated the matrix elements of the four-quark operators relevant to 
the B s meson lifetime difference in QCD sum rules in HQET. The sum rules are constructed 
with three-point correlators and both the perturbative and nonperturbative contribution are 
taken into account. Our result shows that the usual factorization assumption is indeed a 
good approximation. The numerical results show that the sum rules of those operators have 
a good platform. The perturbative contribution to sum rules are about 45 — 65% of the total 
contribution. Our results are not sensitive torn,. The life difference „ Ba is found to be around 

(7.0 ± 0.8) x 10~ 2 . This result is compatible with those predicted by lattice calculations. The 
a s corrections have not been taken into account in the sum rules and they will definitely have 
effects on the resulting numerical values. To get more accurate prediction, the a s corrections 
should be taken into account, which is beyond the content of the letter. 
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Fig. 1 



Figure 1: Dominant non-vanishing Feynman diagrams for T°(u,u') 




Fig2 
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Figure 2: The dependence of (B s \Ol\B s ) on T,u c 
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u, = 0.S5 GeV 
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jure 3: The variation of \(B s \Os\B s )\ with T, oj c 




Figure 4: The variation of 1Z with T, oo c 
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